On Pattern-Cladistic Analyses Based on Complete Plastid Genome Sequences.

The fundamental Hennigian principle, grouping solely on synapomorphy, is seldom used in modern phylogenetics. In the submitted paper, we apply this principle in reanalyzing five datasets comprising 197 complete plastid genomes (plastomes). We focused on the latter because plastome-based DNA sequence data gained dramatic popularity in molecular systematics during the last decade. We show that pattern-cladistic analyses based on complete plastid genome sequences can successfully resolve affinities between plant taxa, simultaneously simplifying both the genomic and analytical frameworks of phylogenetic studies. We developed "Matrix to Newick" (M2N), a program to represent the standard molecular alignment of plastid genomes in the form of trees or relationships directly. Thus, massive plastome-based DNA sequence data can be successfully represented in a relational form rather than as a standard molecular alignment. Application of methods of median supertree construction (the Average Consensus method has been used as an example in this study) or Maximum Parsimony analysis to relational representations of plastome sequence data may help systematist to avoid the complicated assumption-based frameworks of Maximum Likelihood or Bayesian phylogenetics that are most used today in massive plastid sequence data analyses. We also found that significant amounts of pure genomic information that typically accommodate the majority of current plastid phylogenomic studies can be effectively dropped by systematists if they focus on the pattern-cladistics or relational analyses of plastome-based molecular data. The proposed pattern-cladistic approach is a powerful and straightforward heuristic alternative to modern plastome-based phylogenetics.

Akhania, a new genus for Salsoladaghestanica, Caroxyloncanescens and C.carpathum (Salsoloideae, Chenopodiaceae, Amaranthaceae).

Genus Salsola s.l. was recently split into several genera of different phylogenetic placements within Salsoloideae, but both taxonomic and phylogenetic relationships of some parts of the former broadly defined Salsola still need to be clarified. A remarkable example is Salsolacanescens nom. illegit. ≡ Salsolaboissieri, a taxon with tricky taxonomic history that was only recently transferred to the genus Caroxylon (tribe Caroxyleae). Salsoladaghestanica, a narrow endemic of Central Dagestan (Russian Federation), was not even included in previous molecular studies of Salsoloideae and therefore still lacks an appropriate estimation of its relationships. Molecular phylogeny constructed here using nuclear and plastid DNA sequence data clearly placed Salsoladaghestanica and Caroxyloncarpathum as sister taxa and the clade S.daghestanica, Caroxyloncanescens (Salsolaboissieri), C.carpathum (Salsolacarpatha) as a sister of the monophyletic Caroxylon. All three species are distinct from Caroxylon from a morphological standpoint. In conclusion, a new genus, Akhania, was established for these taxa. The detailed distribution of Akhaniadaghestanica is presented for the first time.

A new, simple, highly scalable, and efficient protocol for genomic DNA extraction from diverse plant taxa.

PREMISE: Commonly used molecular techniques such as next-generation sequencing require reliable methods to extract DNA quickly and efficiently. Secondary compounds within plant tissues make this requirement all the more challenging, often forcing researchers to test different extraction methods tailored to their particular species of interest in order to obtain large amounts of high-quality genomic DNA. The opportunities provided by high-throughput, next-generation sequencing only exacerbate these problems, especially when trying to extract DNA from multiple species at the same time. Several methods have attempted to resolve the challenges of obtaining suitable DNA from plants; however, a rapid, high-yield, high-quality, and highly scalable DNA extraction method is still needed. METHODS AND RESULTS: We present a rapid DNA extraction protocol that utilizes a buffer with relatively large amounts of cetyltrimethylammonium bromide (CTAB) and sodium chloride, combined with a silica maxi-column cleanup of the extracted DNA. The new method is easy to implement using standard equipment and inexpensive reagents. The entire procedure (from grinding to the final elution) can be completed in less than two hours for a single sample and can be easily scaled to meet desired research goals. It works on diverse green plants with highly varied secondary chemistries (e.g., ferns, gymnosperms, and phylogenetically divergent angiosperms). CONCLUSIONS: Application of the protocol to various plant species yielded DNA of high quality in less than two hours and can be adjusted to extract DNA at large (maxi-preps) or small (96-well minipreps) scales. We anticipate that our method will be of wide utility for rapidly isolating large quantities of quality genomic DNA from diverse plant species and will have broad applications in phylogenetic studies utilizing PCR and short-read DNA sequencing.

Synapomorphies Behind Shared Derived Characters: Examples from the Great Apes' Genomic Data.

Phylogenetic systematics (e.g., cladistics) is one of the most important analytical frameworks of modern Biology. It seems to be common knowledge that within phylogenetics, 'groups' must be defined based solely on the synapomorphies or on the "derived" characters that unite two or more taxa in a clade or monophyletic group. Thus, the idea of synapomorphy seems to be of fundamental influence and importance. Here I will show that the most common and straightforward understanding of synapomorphy as a shared derived character is not sufficient and eventually must be rejected in favor of Nelson's relational interpretation of such term. Arguing for this point and using three examples from previously published Apes' genomic matrices, I explicitly demonstrate that the relationship (Pongo (Gorilla (Homo, Pan))) with Hylobatidae as a sister taxon, may be successfully recovered by three-taxon statement analysis (3TA) and three-taxon statement average consensus analysis (3TS-ACA) even if all of the evident standard shared derived molecular characters of the relationship (Pongo (Gorilla (Homo, Pan))) with Hylobatidae as a sister taxon, have been excluded from the molecular alignments. Neither conventional Maximum Parsimony nor Maximum Likelihood or Bayesian Inference can do this in such situation. Thus, our results show that the relationship (Pongo (Gorilla (Homo, Pan))) with Hylobatidae as a sister taxon has appeared, in some way, behind standard shared derived characters: the last ones could be excluded, but the relationship remains the same.

Temporal area approach for distributional data in biogeography.

A structural approach to temporality in distributional data for use in palaeobiogeography is described herein. Pre-established areas in the distributional data matrix are split temporally, allowing a single geographical space to have multiple iterations [e.g. Area A (Lower Devonian), Area A (Middle Devonian)]. The resulting temporal matrix will allow the representation and capture of any differing relationships through time. Designed primarily for Parsimony Analysis of Endemicity (PAE) and biotic similarity analyses, this approach simply structures distributional data within a temporal partition, meaning that numerical methods can be used to assess relationships between areas to find a branching diagram. Created through the application of the temporal matrix to a given analysis, Temporal Area Approach (TAAp) is a structural approach that facilitates exploration of the data rather than being a hypothesis-driven model following analysis. Understanding the behaviour of non-phylogenetic palaeobiogeographical data and reducing the prevalence of temporal artefacts will lead to more robust area classifications.

Application of CRISPR/Cas9 to Tragopogon (Asteraceae), an evolutionary model for the study of polyploidy.

Tragopogon (Asteraceae) is an excellent natural system for studies of recent polyploidy. Development of an efficient CRISPR/Cas9-based genome editing platform in Tragopogon will facilitate novel studies of the genetic consequences of polyploidy. Here, we report our initial results of developing CRISPR/Cas9 in Tragopogon. We have established a feasible tissue culture and transformation protocol for Tragopogon. Through protoplast transient assays, use of the TragCRISPR system (i.e. the CRISPR/Cas9 system adapted for Tragopogon) was capable of introducing site-specific mutations in Tragopogon protoplasts. Agrobacterium-mediated transformation with Cas9-sgRNA constructs targeting the phytoene desaturase gene (TraPDS) was implemented in this model polyploid system. Sequencing of PCR amplicons from the target regions indicated simultaneous mutations of two alleles and four alleles of TraPDS in albino shoots from Tragopogon porrifolius (2x) and Tragopogon mirus (4x), respectively. The average proportions of successfully transformed calli with the albino phenotype were 87% and 78% in the diploid and polyploid, respectively. This appears to be the first demonstration of CRISPR/Cas9-based genome editing in any naturally formed neopolyploid system. Although a more efficient tissue culture system should be developed in Tragopogon, application of a robust CRISPR/Cas9 system will permit unique studies of biased fractionation, the gene-balance hypothesis and cytonuclear interactions in polyploids. In addition, the CRISPR/Cas9 platform enables investigations of those genes involved in phenotypic changes in polyploids and will also facilitate novel functional biology studies in Asteraceae. Our workflow provides a guide for applying CRISPR/Cas9 to other nongenetic model plant systems.

10KP: A phylodiverse genome sequencing plan.

Understanding plant evolution and diversity in a phylogenomic context is an enormous challenge due, in part, to limited availability of genome-scale data across phylodiverse species. The 10KP (10,000 Plants) Genome Sequencing Project will sequence and characterize representative genomes from every major clade of embryophytes, green algae, and protists (excluding fungi) within the next 5 years. By implementing and continuously improving leading-edge sequencing technologies and bioinformatics tools, 10KP will catalogue the genome content of plant and protist diversity and make these data freely available as an enduring foundation for future scientific discoveries and applications. 10KP is structured as an international consortium, open to the global community, including botanical gardens, plant research institutes, universities, and private industry. Our immediate goal is to establish a policy framework for this endeavor, the principles of which are outlined here.

A laid-back trip through the Hennigian Forests.

BACKGROUND: This paper is a comment on the idea of matrix-free Cladistics. Demonstration of this idea's efficiency is a major goal of the study. Within the proposed framework, the ordinary (phenetic) matrix is necessary only as "source" of Hennigian trees, not as a primary subject of the analysis. Switching from the matrix-based thinking to the matrix-free Cladistic approach clearly reveals that optimizations of the character-state changes are related not to the real processes, but to the form of the data representation. METHODS: We focused our study on the binary data. We wrote the simple ruby-based script FORESTER version 1.0 that helps represent a binary matrix as an array of the rooted trees (as a "Hennigian forest"). The binary representations of the genomic (DNA) data have been made by script 1001. The Average Consensus method as well as the standard Maximum Parsimony (MP) approach has been used to analyze the data. PRINCIPLE FINDINGS: The binary matrix may be easily re-written as a set of rooted trees (maximal relationships). The latter might be analyzed by the Average Consensus method. Paradoxically, this method, if applied to the Hennigian forests, in principle can help to identify clades despite the absence of the direct evidence from the primary data. Our approach may handle the clock- or non clock-like matrices, as well as the hypothetical, molecular or morphological data. DISCUSSION: Our proposal clearly differs from the numerous phenetic alignment-free techniques of the construction of the phylogenetic trees. Dealing with the relations, not with the actual "data" also distinguishes our approach from all optimization-based methods, if the optimization is defined as a way to reconstruct the sequences of the character-state changes on a tree, either the standard alignment-based techniques or the "direct" alignment-free procedure. We are not viewing our recent framework as an alternative to the three-taxon statement analysis (3TA), but there are two major differences between our recent proposal and the 3TA, as originally designed and implemented: (1) the 3TA deals with the three-taxon statements or minimal relationships. According to the logic of 3TA, the set of the minimal trees must be established as a binary matrix and used as an input for the parsimony program. In this paper, we operate directly with maximal relationships written just as trees, not as binary matrices, while also using the Average Consensus method instead of the MP analysis. The solely 'reversal'-based groups can always be found by our method without the separate scoring of the putative reversals before analyses.

A broad-scale comparison of aerobic activity levels in vertebrates: endotherms versus ectotherms.

Differences in the limits and range of aerobic activity levels between endotherms and ectotherms remain poorly understood, though such differences help explain basic differences in species' lifestyles (e.g. movement patterns, feeding modes, and interaction rates). We compare the limits and range of aerobic activity in endotherms (birds and mammals) and ectotherms (fishes, reptiles, and amphibians) by evaluating the body mass-dependence of VO2 max, aerobic scope, and heart mass in a phylogenetic context based on a newly constructed vertebrate supertree. Contrary to previous work, results show no significant differences in the body mass scaling of minimum and maximum oxygen consumption rates with body mass within endotherms or ectotherms. For a given body mass, resting rates and maximum rates were 24-fold and 30-fold lower, respectively, in ectotherms than endotherms. Factorial aerobic scope ranged from five to eight in both groups, with scope in endotherms showing a modest body mass-dependence. Finally, maximum consumption rates and aerobic scope were positively correlated with residual heart mass. Together, these results quantify similarities and differences in the potential for aerobic activity among ectotherms and endotherms from diverse environments. They provide insights into the models and mechanisms that may underlie the body mass-dependence of oxygen consumption.

What is Atraphaxis L. (Polygonaceae, Polygoneae): cryptic taxa and resolved taxonomic complexity instead of the formal lumping and the lack of morphological synapomorphies.

BACKGROUND: The recently proposed recircumscription of the genus Atraphaxis (incl. Atraphaxis section Ovczinnikovia O.V. Yurtseva ex. S. Tavakkoli and Polygonum sect. Spinescentia Boissier (=A. sect. Polygonoides S. Tavakkoli, Kaz. Osaloo & Mozaff.) makes this genus fairly heterogeneous and therefore almost undefinable based on morphology. A critical comprehensive reappraisal of the group is necessary. METHODS: Using the DNA sequence data (ITS1&2 regions of nrDNA and combined trnL intron + trnL-F IGS and rpl32-trnL((UAG)) IGS regions of plastid genome), Maximum Likelihood (ML) and Bayesian analyses (BI) were applied for phylogenetic reconstructions of the tribe Polygoneae with special attention to Atraphaxis, and related taxa. Maximum parsimony reconstructions of the evolution of perianth morphology and sporoderm ornamentation in the tribe Polygoneae were also performed. Life history, morphology of shoots, leaf blades, ocreas, perianth and achene morphology, ultrasculpture of achene surface, and pollen morphology were compared, and SEM and LM images were provided. PRINCIPAL FINDINGS: The genera Atraphaxis and Polygonum were found to be widely polyphyletic. The rarest and morphologically remarkable endemic of Tian-Shan and Pamir Atraphaxis ovczinnikovii (Atraphaxis sect. Ovczinnikovia O.V. Yurtseva ex. S. Tavakkoli) was confirmed to be a sister of the clade (Atraphaxis + Polygonum sect. Spinescentia) in plastid topology. The genus Bactria (=Atraphaxis sect. Ovczinnikovia), which circumscribes two species, is newly established as a result of this analyses. Morphological data confirm the originality of the taxon. DISCUSSION: We are arguing for a narrow delimitation of Atraphaxis with petalloid segments and striato-perforate sporoderm ornamentation as morphological synapomorphies. The recently proposed inclusion of Polygonum sect. Spinescentia in Atraphaxis is fairly questionable from a morphological standpoint. The rank of Polygonum sect. Spinescentia requires further clarification. The generic composition of the tribe Polygoneae also requires future reappraisals.

Body mass scaling of passive oxygen diffusion in endotherms and ectotherms.

The area and thickness of respiratory surfaces, and the constraints they impose on passive oxygen diffusion, have been linked to differences in oxygen consumption rates and/or aerobic activity levels in vertebrates. However, it remains unclear how respiratory surfaces and associated diffusion rates vary with body mass across vertebrates, particularly in relation to the body mass scaling of oxygen consumption rates. Here we address these issues by first quantifying the body mass dependence of respiratory surface area and respiratory barrier thickness for a diversity of endotherms (birds and mammals) and ectotherms (fishes, amphibians, and reptiles). Based on these findings, we then use Fick's law to predict the body mass scaling of oxygen diffusion for each group. Finally, we compare the predicted body mass dependence of oxygen diffusion to that of oxygen consumption in endotherms and ectotherms. We find that the slopes and intercepts of the relationships describing the body mass dependence of passive oxygen diffusion in these two groups are statistically indistinguishable from those describing the body mass dependence of oxygen consumption. Thus, the area and thickness of respiratory surfaces combine to match oxygen diffusion capacity to oxygen consumption rates in both air- and water-breathing vertebrates. In particular, the substantially lower oxygen consumption rates of ectotherms of a given body mass relative to those of endotherms correspond to differences in oxygen diffusion capacity. These results provide insights into the long-standing effort to understand the structural attributes of organisms that underlie the body mass scaling of oxygen consumption.

250 years of hybridization between two biennial herb species without speciation.

Hybridization between plant species can generate novel morphological diversity and lead to speciation at homoploid or polyploid levels. Hybrids between biennial herbs Tragopogon pratensis and T. porrifolius have been studied in experimental and natural populations for over 250 years. Here we examine their current status in natural populations in southeast England. All hybrids found were diploid; they tended to grow taller and with more buds than their parental species; many showed partial fertility; a few showed evidence of backcrossing. However, we found no evidence to suggest that the hybrids are establishing as a new species, nor can we find literature documenting speciation of these hybrids elsewhere. This lack of speciation despite at least 250 years of hybridization contrasts with the fact that both parental species have formed new allopolyploid species through hybridization with another diploid, T. dubius. Understanding why hybrids often do not speciate, despite repeated opportunities, would enhance our understanding of both the evolutionary process and risk assessments of invasive species.

One-seeded fruits in the core Caryophyllales: their origin and structural diversity.

The core Caryophyllales consist of approximately 30 families (12,000 species) distributed worldwide. Many members evolved one-seeded or conjoined fruits, but their origin and structural diversity have not been investigated. A comparative anatomical investigation of the one-seeded fruits within the core Caryophyllales was conducted. The origin of the one-seeded fruits and the evolutionary reconstructions of some carpological characters were traced using a tree based on rbcl and matK data in order to understand the ancestral characters and their changes. The one-seeded fruit type is inferred to be an ancestral character state in core Caryophyllales, with a subsequent increase in the seed number seen in all major clades. Most representatives of the 'Earlier Diverging' clade are distinguished in various carpological traits. The organization of the pericarp is diverse in many groups, although fruits with a dry, many-layered pericarp, consisting of sclerenchyma as outer layers and a thin-walled parenchyma below, with seeds occupying a vertical embryo position, are likely ancestral character states in the core Caryophyllales clade. Several carpological peculiarities in fruit and seed structure were discovered in obligate one-seeded Achatocarpaceae, Chenopodiaceae, Nyctaginaceae, Seguieriaceae and Sarcobataceae. The horizontal embryo evolved in only certain groups of Chenopodiaceae. The bar-thickening of endotegmen cells appears to be an additional character typical of core Caryophyllales. The syncarpy-to-lysicarpy paradigm in Caryophyllaceae needs to be reinterpreted.

Multiple origins and chromosomal novelty in the allotetraploid Tragopogon castellanus (Asteraceae).

Tragopogon includes two classic examples of recently formed allopolyploid species in North America: T. mirus and T. miscellus. Older Tragopogon allotetraploids from Eurasia offer ideal taxa for comparing the longer term outcomes of allopolyploidy. To help resolve the ancestry of one of these older polyploids, phylogenetic analyses of multiple populations of the allotetraploid T. castellanus (2n = 24) and its putative diploid parents, T. crocifolius and T. lamottei, were conducted using sequences from nuclear (internal transcribed spacer, ITS; and alcohol dehydrogenase 1A, Adh) and plastid (trnT-trnL spacer, trnL intron, trnL-trnF spacer and rpl16 intron) loci. Fluorescence in situ hybridization (FISH) and genomic in situ hybridization (GISH) were used to investigate the chromosomal constitution of T. castellanus. Our data confirm that the widely distributed T. crocifolius and the Iberian endemic, T. lamottei, are the diploid parents of T. castellanus, and that this polyploid formed at least three times. One group of populations of T. castellanus is distinct in exhibiting two pairs of rearranged chromosomes. These data suggest that some of the chromosomal variants that originate in young polyploids (here, an intergenomic translocation) may become fixed in populations, contributing to novelty in older polyploid lineages. The geographical distributions of the allopolyploids and parents are also complex, with allotetraploid populations being disjunct from one or both of the most closely related diploid parental populations.

Phylogeny, divergence times, and historical biogeography of the angiosperm family Saxifragaceae.

Saxifragaceae (Saxifragales) contain approximately 640 species and 33 genera, about half of which are monotypic. Due to factors such as morphological stasis, convergent morphological evolution, and disjunct distributions, relationships within Saxifragaceae have historically been troublesome. The family occurs primarily in mountainous regions of the Northern Hemisphere, with the highest generic and species diversity in western North America, but disjunct taxa are known from southern South America. Here, we integrate broad gene (56 loci) and taxon (223 species) sampling strategies, both the most comprehensive to date within Saxifragaceae, with fossil calibrations and geographical distribution data to address relationships, divergence times, and historical biogeography among major lineages of Saxifragaceae. Two previously recognized main clades, the heucheroids (eight groups+Saniculiphyllum) and saxifragoids (Saxifraga s.s.), were re-affirmed by our phylogenetic analyses. Relationships among the eight heucheroid groups, as well as the phylogenetic position of Saniculiphyllum within the heucheroids, were resolved with mostly high support. Divergence time estimates indicate that Saxifragaceae began to diversify ca. 38.37 million years ago (Mya; 95% HPD=30.99-46.11Mya) in the Mid-Late Eocene, and that the two major lineages, the heucheroids and saxifragoids, began to diversify approximately 30.04Mya (95% HPD=23.87-37.15Mya) and 30.85 Mya (95% HPD=23.47-39.33Mya), respectively. We reconstructed ancestral geographic areas using statistical dispersal-vicariance (S-DIVA). These analyses indicate several radiations within Saxifragaceae: one in eastern Asia and multiple radiations in western North America. Our results also demonstrate that large amounts of sequence data coupled with broad taxon sampling can help resolve clade relationships that have thus far seemed intractable.

At least 23 genera instead of one: the case of Iris L. s.l. (Iridaceae).

BACKGROUND: Iris L. s.l. is one of the most diverse and well-known genera in the Asparagales, with approximately 250-300 circumscribed species and significant economic impact. The taxonomy of the genus has suffered dramatic changes in the last century, particularly in the last decades after the application of molecular techniques. As a result several contrasting systematic arrangements are currently available to taxonomists. Many genera that were split from Iris s.str. in the past, on the basis of morphology (e.g., Hermodactylus, Iridodictyum, Juno, Pardanthopsis, and Xiphion, among others), are now a priori re-included in a very widely circumscribed Iris s.l. (incl. Belamcanda). This resulted in a more heterogeneous genus that is more difficult to define on morphological grounds. Testing congruence between taxonomic treatments and the results of recent molecular studies of Iris has never been performed, mostly due to the lack of proper taxonomic context. RESULTS: We generated several conventional phylogenies for Iris & outgroups using extensive sampling of taxa (187) and characters (10 plastid loci). We demonstrate that the natural history of Iris, written either as conventional molecular phylogenies or, if viewing in the context of the comparative approach, as a nested most parsimonious hierarchy of patterns, appear to be fully congruent with the narrow taxonomical treatment of the genus, restricted to the rhizomatous "bearded" taxa. The resulting topologies place Belamcanda, Pardanthopsis, and Gattenhofia as sisters to Iris s.str. and genus Siphonostylis as sister to Iris s.l. CONCLUSION: The present study clearly justifies the splitting of Iris s.l. into at least 23 genera, 18 of which have already been accepted in the past by numerous authorities. These genera are characterized by unique combinations of partly overlapping morphological characters and biogeography. Moreover, nearly the same entities, which we here recognize at a generic rank, were for centuries frequently referred to by horticulturists as "working-name" groups.

Phylogenetic relationships and character evolution analysis of Saxifragales using a supermatrix approach.

PREMISE OF THE STUDY: We sought novel evolutionary insights for the highly diverse Saxifragales by constructing a large phylogenetic tree encompassing 36.8% of the species-level biodiversity. • METHODS: We built a phylogenetic tree for 909 species of Saxifragales and used this hypothesis to examine character evolution for annual or perennial habit, woody or herbaceous habit, ovary position, petal number, carpel number, and stamen to petal ratio. We employed likelihood approaches to investigate the effect of habit and life history on speciation and extinction within this clade. • KEY RESULTS: Two major shifts occurred from a woody ancestor to the herbaceous habit, with multiple secondary changes from herbaceous to woody. Transitions among superior, subinferior, and inferior ovaries appear equiprobable. A major increase in petal number is correlated with a large increase in carpel number; these increases have co-occurred multiple times in Crassulaceae. Perennial or woody lineages have higher rates of speciation than annual or herbaceous ones, but higher probabilities of extinction offset these differences. Hence, net diversification rates are highest for annual, herbaceous lineages and lowest for woody perennials. The shift from annuality to perenniality in herbaceous taxa is frequent. Conversely, woody perennial lineages to woody annual transitions are infrequent; if they occur, the woody annual state is left immediately. • CONCLUSIONS: The large tree provides new insights into character evolution that are not obvious with smaller trees. Our results indicate that in some cases the evolution of angiosperms might be conditioned by constraints that have been so far overlooked.

TAXODIUM version 1.0: a simple way to generate uniform and fractionally weighted three-item matrices from various kinds of biological data.

An open-access program for generating three-item statement (3TS) matrices from data such as molecular sequences does not currently exist. The recently developed LisBeth package allows for representation of hypotheses of homology among taxa or areas directly as rooted trees or as hierarchies; however, LisBeth is not a standard matrix-based platform. Here we present "TAXODIUM version 1.0" (TAXODIUM), a program designed for building 3TS-matrices from binary, additive (ordered) and non-additive (unordered) multistate characters, with both uniform and fractional weighting of the statements. TAXODIUM also facilitates, for the first time, use of Maximum Likelihood analyses with 3TS matrices, but future implementation of the 3TS analysis in a statistical framework will require more exploration.

A maximum likelihood approach to generate hypotheses on the evolution and historical biogeography in the Lower Volga Valley regions (southwest Russia).

The evolution of the diverse flora in the Lower Volga Valley (LVV) (southwest Russia) is complex due to the composite geomorphology and tectonic history of the Caspian Sea and adjacent areas. In the absence of phylogenetic studies and temporal information, we implemented a maximum likelihood (ML) approach and stochastic character mapping reconstruction aiming at recovering historical signals from species occurrence data. A taxon-area matrix of 13 floristic areas and 1018 extant species was constructed and analyzed with RAxML and Mesquite. Additionally, we simulated scenarios with numbers of hypothetical extinct taxa from an unknown palaeoflora that occupied the areas before the dramatic transgression and regression events that have occurred from the Pleistocene to the present day. The flora occurring strictly along the river valley and delta appear to be younger than that of adjacent steppes and desert-like regions, regardless of the chronology of transgression and regression events that led to the geomorphological formation of the LVV. This result is also supported when hypothetical extinct taxa are included in the analyses. The history of each species was inferred by using a stochastic character mapping reconstruction method as implemented in Mesquite. Individual histories appear to be independent from one another and have been shaped by repeated dispersal and extinction events. These reconstructions provide testable hypotheses for more in-depth investigations of their population structure and dynamics.